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ABSTRACT 

Long-chain hydrocarbon anions C„H~ {n = 4,6,8) have recently been found to be abundant in 
a variety of interstellar clouds. In order to explain their large abundances in the denser (prestel- 
lar/protostellar) environments, new chemical models are constructed that include gas-grain inter- 
actions. Models including accretion of gas-phase species onto dust grains and cosmic-ray-induced 
desorption of atoms are able to reproduce the observed anion-to-neutral ratios, as well as the abso- 
lute abundances of anionic and neutral carbon chains, with a reasonable degree of accuracy. Due 
to their destructive effects, the depletion of oxygen atoms onto dust results in substantially greater 
polyyne and anion abundances in high-density gas (with njja ^ 10^ cm~^). The large abundances 
of carbon-chain-bearing species observed in the envelopes of protostars such as L1527 can thus be 
explained without the need for warm carbon-chain chemistry. The CgH" anion-to-neutral ratio is 
found to be most sensitive to the atomic O and H abundances and the electron density. Therefore, as 
a core evolves, falling atomic abundances and rising electron densities are found to result in increasing 
anion-to-neutral ratios. Inclusion of cosmic-ray desorption of atoms in high-density models delays 
freeze-out, which results in a more temporally-stable anion-to-neutral ratio, in better agreement with 
observations. Our models include reactions between oxygen atoms and carbon-chain anions to pro- 
duce carbon-chain-oxide species CeO, C7O, HCeO and HC7O, the abundances of which depend on 
the assumed branching ratios for associative electron detachment. 
Subject headings: astrochcmistry - ISM: abundances - ISM: clouds - ISM; molecules 



1. INTRODUCTION 

Hydrocarbon anions have been shown to be an im- 
portant part of the molecular inventory of the Galaxy 
by their recent discoveries in the qui escent molecular 
clouds TMC-1 (iMcCarthv et al.l |2Q06| ) and Lupus-lA 
(jSakai et al.ll2010[ ). the prestellar cores L1544 and L1512 
and the envelopes of t he low-mas s prot o stars L1527, 
L1251A and L1521 F (ISakai et al.l 120071: IGupta et al.l 
[2009t ICordiner et all 120111) . Previous chemical models 
have been able to reproduce with reasonable accuracy the 
observed abunda nces of CsH' and CgH' in TMC-1 and 
IRC-H0216 fe o. IMillar et al.l[2007t IRemiian et al.l [2007l: 
ICordiner eFall 120081 : iWalsh et al.ll2009D. but not so well 
in th e denser environment of L1527 (iHarada fc Herbstj 
12008'). The abundance s of C N" and C.^j N- observed in 
IRC+10216 (Thaddeus 'et^^OS : Agu ndez et al.l[20To[ ) 
closely match the predictions m ade by the chemical 
model of ICordiner k Millaii (|2009f ). These studies show 
that the basic gas-phase chemical processes involving 
molecular anions are fairly well understood. So far, 
however, interactions between negatively-charged gas- 
phase species and dust grain surfaces have been largely 
neglected in chemical models for anions. iFlower et al.l 
([2007) modelled the CgH" to CeH abundance ratio 
([C6H^]/[C6H]) in 'typical' dark-cloud gas (with den- 
sity 10"* cm~'^), using a chemical model accounting for 
the detailed calculation of charge balance, including 
gas, dust grain and PAH charge exchange, but did not 
consider the effects of depletion/freeze-out of molecules 
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onto the grains. Gas-grain interactions are theorized 
to be of fundamental importance for the abundances 
of gas-phase molecules iii dense interstellar clou d s {e.g . 
IBrown fc CharnlevI IT991 . and ICordiner et aH mm 
proposed that oxygen freeze-out could help to explain 
the large abundances of carbon-chain-bearing species and 
anions observed in low-mass protostellar envelopes and 
prestellar cores with densities riHa ~ 10^ — 10^ cm"-^. 

The total abundance of molecular anions is believed to 
be an important factor affecting the ionization balance 
of the interstellar mediu m (ISM) (e.g. Lepp & Dalgarn3 
[1981 IFlower et al.llMTt Wakelam & Hcrbst 2008). An- 
ion abundances are theorized to be sensitive to elec- 
tron attachment and photodetachment rates, and may 
therefore provide a useful tool for the determination 
of electron densities and cosmic ray /photo- ionization 
rates in as trophysical environments (jMillar et al.l 1200*71 : 
IFlower et~a l. 2007). As shall become apparent from the 
present study, the interstellar CgH" anion-to-neutral ra- 
tio can also be a tracer of the abundances of gas-phase 
atomic hydrogen, oxygen and carbon, due to the rapid 
reactions between hydrocarbon anions and these species. 

This article presents new chemical models for inter- 
stellar carbon chains and their anions, with the aim of 
understanding the more recent CgH" detections in rel- 
atively dense prestellar cores and cold protostellar en- 
velopes. The inclusion of gas-grain interactions at vary- 
ing levels of complexity permits analysis of the impor- 
tance of these processes for anion chemistry, particularly 
at higher densities. We also discuss the possibility that 
reactions between hydrocarbon anions and oxygen atoms 
can produce detectable abundances of the carbon chain 
oxides C„0 and HC„0 [n = 6, 7) and the smaller hydro- 
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carbon anion C2H ; these species are yet to be detected 
in space. 

2. CHEMICAL MODELLING 

The chemical models presen ted here utilize the (dipole- 
enhanced) RATE06 network (iWoodall et al 112007.) . aug - 
mented with the anion chemistry of IWalsh et alj (|2009l ) , 
and extended further to include CaH" and C2H~. An- 
ions in the model are formed by radiatiye electron at- 
tachment, with ra t es tak en from the calculations by 
iHerbst fc Osamural ()2008[ ). Anions are destroyed in re- 
actions with H, C, N and O atoms and in mutual- 
neutralisation reactions with the most abundant cations 
(including C+, H+, H+ and HCO+). Product branch- 
ing ratios for ani o n-neu tral reactions are taken from 
lEichclbcrg er et all ()2007[ ). with additional data from V. 
Bierbaum (priyate communication, 2010). For reactions 
inyolving CaH", the average was taken of the published 
rate coefficients for C2H~ and C4H^. In addition to 
the C3O molecule already present in the RATE06 net- 
work, we haye added the longer carbon-chain oxides 
C„0, HC„0 (for n = 6 and 7), as well as their pro- 
tonated forms ((H)C„OH+) and anions ((H) C„0^) (see 
Section [5] for details) . PAHs haye not been included in 
our models due to a lack of eyidenc e concerning t heir 
abundances in dark molecular clouds (jPeeters II2011[ ). 

The abundances of 454 gas-phase species are calcu- 
lated as a function of time, linked by 5132 binary reac- 
tions. Models begin at time t = with neutral atomic 
abundances giyen in Table [TJ The abundances of heavier 
metals (Na, Mg, Fe, P etc.) are set to zero, cons istent 
with their non- detection in the ISM (jTurnerill99ll ). The 
cosmic ray ionisation rate is set at 1.3 x 10~^^ and 
the visual extinction (Ay) is set to 10 mag. Hydrogen is 
predominantly molecular and models are run at densities 
= 10^, 10^ and 10^ cm^'^. These span the range of 
interstellar cloud densities where anions have so far been 
observed, from quiescent, dense molecular clouds {e.g. 
TMC-1, with TiHa = 10^ cm~^), to prestellar cores {e.g. 
L1512, with riHa — 10^ cm~'^) and protostellar envelopes 
{e.g. L1527, with riH2 = 10^ cm~^). The kinetic temper- 
ature is fixed at 10 K. These simplified models are useful 
for analysing basic astrochemical processes, and we do 
not consider the impact on the chemistry of protostellar 
heating, nor dynamical effects such as infall or outflows 
that may be present in some of the sources. Throughout 
this article, fractional abundances are given with respect 
to the total H2 number density (rifja)- 

2.1. Model l: No accretion 

Our first model considers only gas-phase chemistry and 
neglects the gas-grain interaction apart from the forma- 
tion of H2 on dust grains, which is assumed to occur 
at a rate 4.7 x 10~^*nH(?^H + 2nH->)- This modelling 
approach is the sam e as that used by IHarada fc Herbstl 
p08) for L1527 and IWalsh et"all ([20091) for TMC-1 and 
is therefore useful as a benchmark for comparison with 
those models. For de tails of the modelUng pr ocedure see 
IWoodall et aH (l2007l) and IWalsh etHI (pOOl) . 

2.2. Model II; Accretion without desorption 

This model builds on Model i with the inclusion of gas- 
grain interactions. These are treated using the approach 



of iCharnlevI ()1997D . whereby gas-phase species (except 
for the highly volatile H+, H2, He and Hj^) that collide 
with dust grains have a certain probability of sticking to 
the surface. We model a dust grain population composed 
of three distinct charge states: G° (neutral), (single 
positive charge) and G~ (single negative charge). Due to 
electrostatic repulsion, anions are assumed not to stick to 
negatively-charged grains. Anions undergo mutual neu- 
tralisation upon collision with positively-charged grains, 
resulting in transfer of an electron from the anion (X^ ) 
to the grain, and returning the neutral species (X) to the 
gas. Positive ions (X"*") undergo mutual- neutralisation in 
collisions with negatively-charged grains, again return- 
ing the neutral species (X) to the gas. The dust grain 
charges are calculated as a function of time using the for- 
mulation of Umcbavashi & Nakano (1990) for a uniform- 
size grain population, with the addition of the anion- 
grain processes mentioned above. The parameters used 
are: dust grain radius = 10^^ cm, dust grain abundance 
— 2.0 X 10~^^, electron sticking-coefficient = 0.6, gas 
sticking-coefficient — 1, temperature T = T^ust — Tgas- 
The charge balance for this chemical plasma is 

n(X-) + n{G-) +ne=J2 "(X^) + ^HG+) (1) 
X- x+ 

where the sums are over the number densities of all 
anions (X~) and cations (X+) in the model. 

Due to the rapidity of collisions of electrons with the 
grains (as compared to cations with grains), the major- 
ity of the dust attains a negative charge soon after the 
start of the model, resulting in an approximately steady 
fraction of ~ 90% of the total number of grains as G~ 
during the chemically-relevant timescales. At a density 
of 10^ cm""^, gas-grain collisions result in near-complete 
depletion of the gas after ~ 10^ yr (see Figure [31 middle 
panel) . 

2.3. Model ill: Accretion with desorption of atoms 

Gas-phase species deposited on the dust grains build 
up over time to form an icy mantle (see for example 
lOberg et al]|2011[ ). Energetic processes are theorized to 
result in t he liberation of m antle species back into the 
gas phase ([Leger et al.lll985D . Thermal desorption is very 
slow at dust temperatures of 10 K, and we assume that 
photo-desorption may be neglected in the dark clouds we 
are considering. However, heating of the dust grains by 
cosmic ray impact results in significantly elevated tem- 
peratures for short time periods, s ufficient for volatile 
speci es to be quite rapidly desorbed ([Hasegawa fc Herbs^ 
|1993() . Earlier models have proposed that atoms may ei- 
ther stick inefficiently on dust grai ns, or are selectively 
desorbed with respect to molecules (iHasega w a fc Herbsti 
fl993tlCharnlev et al.ll200ll:IFlower et al.ll2005l) . In our fi- 
nal model, we utilize the rate coefficients for cosmic-ray- 
in duced desorption of C, N, O and S atoms calculated 
by IHasegawa fc Herbs^ (|1993l ) . For simplicity, and due 
to their expected higher binding energies, desorption of 
molecular species is not considered. 

3. RESULTS 

For each of the three models, fractional abundances 
of the carbon-chain-bearing species HC3N, C4H, CeH 
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TABLE 1 

Initial fractional 
abundances 



Species 


Abundance^ 


H 


1.0x10-3 


He 


2.8x10-1 


C 


1.2x10-4 


N 


4.0x10-5 





3.2x10-* 


s 


2.0x10-** 



^ Relative to njja ■ 



and anions C4H and CgH are plotted as a func- 
tion of time in Figures [T][3] at densities riHg = iC, 
10^ and 10^ cm-'^, respectively. Horizontal broken 
lines show the observed fractional abundances of each 
species: for = lO'' cm^'^, TMC-1 data are shown 
and for nB_2 — 10^ cm^'^, L1512 data are shown. For 
riHa — 10^ cm"'^, observed abundances are plotted for 
three sources: the pre-stellar core L1544 (dashed lines), 
the very low-luminosity, young protostar L1521F (dot- 
ted lines), and the Class O/I protostar L1527 (dot-dashed 
lines). We briefly examine the results obtained for each 
density regime. 

3.1. — 10^ cmT^ - Quiescent molecular cloud 

For riHa = 10'* cm"'^ (Figure [J), Models i-iii achieve 
a reasonable match with observations of the quiescent 
molecular cloud TMC-1 (within about an order of mag- 
nitude) at times around 10^ yr, consistent with previ- 
ous models for th is source (e.g. IHerbst fc Lemli 11981 
I Walsh et al]|2009[ ). At times later than this, the HC3N 
and C4H abundances fall well below their respective ob- 
served values. 

Interstellar polyyne (Cn H) and cyanqpolyyn e (HC„N) 
chemistry was discussed bv lWalsh et al.l (|2009D : polyynes 
are radical species and react rapidly with oxygen atoms 
whereas the closed-shell cyanopolyynes do not. Thus, 
in Model 11 the freezing-out of atomic oxygen onto dust 
grains eliminates one of the main polyyne destruction 
channels which significantly extends their lifetime in the 
gas. HC3N, on the other hand, is primarily destroyed by 
reaction with C"*" and continues to fall in a similar man- 
ner between 10^ and 10^ yr independently of whether 
0-atoms are accreting out (see Figure [Tl upper and mid- 
dle panels) . The cosmic-ray-desorption of 0-atoms that 
occurs in Model ill (Figure [1] lower panel) results in sim- 
ilar abundances to Model i because at this relatively low 
density the gas accretion rate is sufficiently low that it is 
counteracted by cosmic-ray desorption, which keeps the 
majority of material off the dust and in the gas-phase. 

At a time of 10^ yr, the modelled CeH- anion-to- 
neutral ratio [CgH-J/fCgH] is equal to 2.3% in all three 
models with rifja — 10^ cm"^ and is thus apparently 
quite insensitive to gas-grain interactions at this time. 
The anion-to-neutral ratio rises at later times, predomi- 
nantly as a result of rising electron densities and falling 
atomic H, C and O abundances (see Section [B]). Com- 
parison of the modelled anion-to-neutral ratio with the 
observed values given in Table [2] shows good agreement 
with the quiescent clouds TMC-1 and Lupus lA. 
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Fig. 1. — Carbon chain and anion fractional abundances as a 
function of time for rtHj = 10* cni~^. Results of Model I (no 
accretion) are shown in the top panel, Model II (accretion with 
no desorption) in the middle panel and Model III (accretion + 
CR-desorption) in the bottom panel. TMC-1 ob served abun- 
dances are shown with dashed lines (data from iSuzuki et al.lll992l : 
IBriinken eFa l. 2007; Agundcz ct al. 200l, with iV(H2) = 1.0 x 10^2 
cm-2 from lCernicharo fc Guelinli 19871) . 

3.2. = 10^ cm~^ - Prestellar core 

At riHa ~ 10^ cm^"^ (Figure[2]), a reasonable agreement 
is obtained with observations of the prestellar core L1512 
at times around (6-8) x 10^ yr. The polyynes and anions 
in Models i and ill (top and bottom panels) show a good 
fit to observations at this time (within a factor of a few) , 
but HC3N is significantly over-predicted. On the other 
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Fig. 2. — As Figure[T]but for nji^ = 10^ cm~^. L1512 observed 
abundances are shown with dashed Hnes (data from'Cordi ner et al.l 
[20T1I. with Af(H2) = 2.0 X 10^2 cm-^ from Kirk ct al. 200D. 

hand, in Model 11 (with freeze-out and no desorption), 
HC3N matches observations at about 6 x 10'' yr whereas 
CgH and CgH" are over-predicted by about an order of 
magnitude; the abundances of the latter are greater as a 
result of 0-atom freeze-out. 

The best fit to the observed CgH^ anion-to-neutral 
ratios in the sources with = 10^ cm""^ (L1512 and 
L1251A; see Table[2]) is also obtained in Models i and ill, 
at around 7 X 10"^ yr, at which time [CeH-J/peH] = 3.1% 
and 3.5% respectively. The freeze-out model with no des- 
orption (Model 11), significantly over-predicts the anion- 




Tlme (yr) 

Fig. 3. — As Figure [1] but for riHa = 10® cm~^. L1544 abun- 
dances are shown with dashed lines (data from [Suzuki et aTll 19921 
and IGupta et al.l pOOS . with N(H2) = 9.4 X lO^^ cm-^ from 
ICrapsi et al.l 12003). L152 1F abundances are shown with dot- 
dashed hnes l|Gupta et aLI [2)09, and our own unpubhshed On- 
sala 3 mm HC3N data - for details of the Onsala observations 
sec Buckle ct al. 2006. We used Af(H2) = 1.4 X 10^^ cm'^ from 
[Crapsi ct al. 2005). L1527 abundances plotted with dotted lines 
(j0rgcnscn ct al. 2004; Sakai ct al. 2007; Agundcz ct al. 2008, with 
Ar(H2) = 9.4 X lO^^ cm"^ f^.^^ Jargcn sen et al.ii2002i ). The L1527 
C4H^ abundance (cyan dotted line) is 2.7 X 10~^^ and lies just 
above the L1521F CsH^ abundance (green dot-dashed line). 
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TABLE 2 

Observed anion-to-neutral ratios, densities and CO 

DEPLETION FACTORS IN DENSE CLOUD CORES 



Source [CeH-J/peH] 


(cm-3) 5(CO) 


Refs. 


TMC-1 


1.6 ±0.3% 


~ 10* 


1.2=^ 


1, 2, 3 


Lupus- 1 A 


2.1 ±0.6% 






4 


L1512 


4.2 ± 1.4% 


~ 10^ 


3.9!^ 


5, 6, 7 


L1251A 


3.6 ± 1.3% 


~ 10^ 




5, 8 


L1544 


2.5 ± 0.8% 


~ 10^ 


7-14 


9, 10, 11 


L1521F 


4.0 ± 1.0% 


~ 10** 


12-15 


9, 10, 12 


L1527 


9.3 ± 2.9% 


~ lO** 


2 4a 


13, 14 


References 




[Brunken et al.' 


(1200' 


3); (2) - 


[Hirahara ct al. 


(1991); (3) 


- Ohishi ct al. 


(199^ 


J); (4) - 



feakai et al., 12010 '): (5) -ICordiner ct al. (2011); (6) - Kirk et al.l 
II2005I ) ; (7) - |. Lcc ct al. (2003); (8) - Lcc ct al. (2010); 
(9) - IGupta et al. (2009,); (10) - Crapsi ct al. (2005); (11) 
- Cascllictal. (200^); (12) - Crapsi ct al. (2004); (13) - 

[Sakai ct al. (2007); (14) - |j0rgen scn ct al. (2004). 



^ Calculated from CO observations assuming an undcpletcd 
abundance of 9.5 X lO^^nua- 

The ILee et"all (2003) CO depiction factor in L1512 has been 
scaled by 0.35 to account for the greater undepleted CO abun- 
dance of 2.7 X lO^^riHj used in that study. 

to-neutral ratio at this time ([C6H~]/[C6H] = 8%) due 
to the larger ratiative electron attachment rate that re- 
sults from the greater electron density that occurs during 
freeze-out (see Section [6l). 

3.3. — 10^ cmT^ - Dense prestellar core / cold 
protostellar envelope 

Gas-grain interactions are more rapid at higher densi- 
ties, so at rtHa = 10^ cm~'^ the effects of accretion on the 
modelled abundances are quite significant, as can be seen 
by comparison of the three panels in Figure[3] Gas-phase 
chemistry also proceeds more rapidly at high densities: 
without accretion (top panel), the conversion of C into 
CO (through the reaction of C and O2), and the presence 
of abundant gas-phase 0-atoms that react with and de- 
stroy the polyynes results in steadily diminishing carbon 
chain abundances over time. Accordingly, the Model i 
abundances fail to match observations of L1544, L1521F 
or L1527 at any time. 

The inclusion of gas-grain interactions results in sub- 
stantially greater carbon chain abundances for t > 
1000 yr in Model 11 compared with Model L Indeed, 
Model II shows a good fit (within an order of magni- 
tude), to the observed abundances in L1544, L1521F and 
L1527 &t t ^ 10^ yr. The only species that does not fit 
well with observations is C4H~ , which has been over- 
predicted by at least an order of mag nitude in all previ- 
ous chemical models for this species CMillar et al .ll2007l: 
Flower et al. 2007; Harada & Herbst 2008; Wals h et al l 
2001 fCordiner fc Millaii 120(391) . The most likely expla- 
nation for this discrepancy is a problem with the theo- 
retical radiative electron attachment rate for C4H, which 
was calculated using a relatively sin iplistic phase-space 
appro ach (for further discussion see iHerbst fc Osamural 

Hoos^. 

The inclusion of cosmic-ray- induced desorption (Model 
III) diminishes the polyyne and anion abundances some- 
what compared with Model 11, primarily due to the in- 
creased gas-phase abundance of destructive 0-atoms. 
The peak that occurs in the polyyne abundances in 
Model II towards the time of freeze-out is eliminated, but 



a good fit to observations is still achieved at ^ 10^ yr, 
although HC3N is again somewhat over-abundant. The 
CgH" anion-to-neutral ratio is found to be quite sensitive 
to the effects of gas-grain interactions at this density; in 
the absence of accretion (Model i), [C6H"]/[C6H] = 1.3% 
at 10* yr (significantly less than observed values; see Ta- 
ble [2|). With accretion and no desorption (Model 11), 
[C6H^]/[C6H] = 8% and rises sharply as time goes on 
(see Figure [5|) . When cosmic-ray-induced desorption of 
atoms is accounted for (Model ill), [C6H^]/[C6H] levels 
off after ~ 10"* yr at about 3%, which matches well with 
the observations of L1544 and L1521F, but somewhat 
less well with the greater value of 9.3% in L1527. 

Dissociative electron attachment (DEA) to the lin- 
ear carbene C6H2 has been suggested as a potentially 
important route to the formation of CgH^ in L1527 
(jSakai et al . 2007). This reaction was shown to be likely 
to proc eed rapidly in the ISM by the theoretical calcula- 
tions of lHerbst fc Osamural (|2008| ). In our models, C6II2 
is included as the more stable poly-acetylene isomer, and 
reactions leading to the carbene (s) are neglected due to 
a lack of laboratory data. The observed (linear carbene) 
C6H2 abundances in TMC-1 and L1527 are reproduced if 
about 1-2% of the total C6H2 in our models is assumed 
to be in the form of the carbene isomer. If the poly- 
acetylene to carbene ratio is fixed at 1%, inclusion of the 
reaction 

CgHs + e- ^ CeH- +H (2) 

with a rate coefficient of 4 .98xl0-'^(r/300)~°-^ 
crn^s"^ (Herb st fc Osamurall2008l ) . generally raises the 
CeH^ abundances in our models by only a negligibly 
small amount. 

C6II2 is destroyed predominantly by reacting with 
atomic carbon to produce C7H -|- H. Its abundance is 
therefore sensitive to the gas-phase C-atom abundance. 
At riHa = 10^ cm~'^ in Model 11 (where C atoms are less 
abundant than in any of the other models), the C6II2 
abundance can become sufficiently large (1.6 x 10~^ at 
7 X 10'' yr) that DEA becomes significant, contributing 
about 50% to the total CeH" formation rate. Observa- 
tions of C6H2 in L1512 and other clouds of similar density 
will be required to further examine the possible signifi- 
cance of DEA for interstellar anion production. DEA is 
not expected to proceed for carbenes smaller than C6II2 
because the reaction becomes endothermic, with an en- 
ergy barrier exceed ing the electron temperat ure in dense 
interstellar clouds (jHerbst fc Osamurall2008l ). 

4. C2H- 

The C2ll~ radiative electron attachment rate is cal- 
culated to be very small (2 x 10-i5(T/300)-°-5 cm^s'^ 
IHerbst fc Osamural[2008[ ) . As such, little C2H" has pre- 
viously been expected in the ISM. However, the dom- 
inant production mechanism in our models is via the 
rapid reaction 

C3H- + O ^ C2H- + CO (3) 

the rate of which is derived from interpolation of 
the C4H~ -I- O a nd C2H^ + O rates measured by 
lEichelberger et al.l (|2007. 1 (and divided by two to account 
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for the 50% branching ratio adopted for the associative 
electron detachment channel). 

In Model ill, the calculated C2H~ fractional abun- 
dance is 6.6xl0~^^ a.t t = 10^ yr for riHa — 10'' cm~'^ 
and 5.9x10"" at t = 10"* yr for nn^ = 10^ cm-^. 
These abundances are low, but may be sufficient to per- 
mit detection with c urrent or future radio telescopes. 
lAgundez et"all (|2008[ ) obtained upper limits on the C2H 
column densities in TMC-1 and L1527 of 2.2 x 10" 
and 1.8 X 10^*^ cm~^ respectively, which correspond to 
fractional abundance upper limits of 2.2 x 10^^^ and 
3.0 X 10^1^ A deeper search for C2H~ in TMC-1 and 
L1527 would thus provide a useful test of our theory with 
regard to anion reactions with 0-atoms. 

5. CARBON CHAIN OXIDES 

In t he laboratory experiments of lEichelberger et aD 
()2007[) on the rate of reaction of atomic oxygen with 
C„(H)~, a loss of ion signal was observed over time (V. 
Bierbaum and Y. Zhibo, private communication, 2010). 
This effect is attributable to associative electron detach- 
ment (AED) reactions 

+ C„(H)- ^ (H)C„0 + e- (4) 

The raw, unpublished laboratory data are consistent 
with a branching ratio for AED of approximately 0.5. 
It must be emphasised, however, that AED products 
were not directly observed, and the existence of the im- 
plied carbon-chain-oxide product channels can only be 
definitively confirmed by additional dedicated labora- 
tory studies. We assume the alternative product channel 
(with braching ratio 0.5) to result in formation of a CO 
molecule and a smaller anion (as described in the previ- 
ous section on C2ll~): 

+ C„(H)- ^ (H)C,-_i + CO (5) 

Model calculations were performed using an AED 
branching ratio of 0.5 in the reaction of oxygen with 
the most abundant anions in the model: C^^, C7, 
CgH- and C7H-, to produce CgO, C7O, HCgO and 
IIC7O, respectively. In addition to the new carbon- 
chain-oxide neutrals, their anions, cations and proto- 
nated forms were also added to the model. Radiative 
electron attachment rates for CeO and C7O were cal- 
culated assuming immediate formation of the stable an- 
ion (through vibronic transitions) , fo howing s-wave elec- 
tron c apture based on the method of lHerbst &: Qsamural 
((2OO8I) . The de generacies of the ground states of 
C„0 for n odd and even are 1 and 3 respectively, 
and 4 for the C„0~ a nions (based on the electroni c 
structure calculations of iRienstra-Kira cofe et aLll2000[ ). 
which leads to radiative electron attachment rates of 
3.26 X 10-'^(T/300)^"-5 cm^s^i for CeO and 9.78 x 
10-'^(T/300)-°-^ cm3 s-i for C7O. In the absence of de- 
tailed electronic structure information for the HCgO" 
and HC70~ anions, an approximate radiative attach- 
ment rate coefficient of 10^''(T/300)~"-^ cm^s"^ has 
been used. 

Based on the reactions and rate coefficients of the 
structurally similar cumulene oxide C3O in the RATE06 
database, the new carbon-chain-oxides have been as- 
sumed to to undergo charge transfer reactions with H''^ 
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Fig. 4. — Modelled carbon-chain-oxide fractional abundances as 
a function of time from Model in at riHj = 10* cm""^. 

and C+, dissociative charge-transfer with He+ and C''' 
and protonation by H3O+ and IICO+. Their cations 
have been assumed to react with II2 and to undergo dis- 
sociative electron recombination. For the recombination 
of HCnO"*" and H2C„0"'", an H atom is assumed to be 
ejected, whereas for CnO"*", a CO molecule is assumed 
to be ejected. Due to the low fractional abundances of 
these ions (~ 10~^^), assuming alternative recombina- 
tion products {e.g. HCnO"*" -I- e~ — > C(„_i)II -I- CO), 
has a negligible effect on the abundances of the species 
of interest in this study. Reactions of (H)C„0^ with H 
and C atoms (modelled using generic rate coefficients of 
10~^ cm-^s"^), are assumed to result in associative de- 
tachment (as ob served in the laboratory f or hydrocarbon 
chain anions bv lEichelberger et al.l 120071 ): their reaction 
with 0-atoms is assumed to result in the production of 
a CO2 molecule and an anion. The (H)C„0~ anion de- 
struction mechanisms are through photodissociation and 
mutual neutralisation with the most abun dant cations 
at a rate 7.51 x 10-^(T/300)-°-5 cm^s"! (|Walsh et all 
l2009f ). Thus, although additional reactions are likely, we 
anticipate that the most rapid processes involving these 
new carbon-chain-oxide species have been - at least ap- 
proximately - accounted for. 

The abundances of the carbon chain oxides C7O, 
IIC7O and their anions (calculated at = 10^ cm^'^ 
using Model III), are plotted as a function of time in 
Figure m The calculated abundances of these modelled 
oxides and their anions are given for densities applica- 
ble to TMC-1 and L1544 in Table El Around the time 
of best agreement of the tih, = 10"* cm^'^ results with 
observations of TMC-1 (- 10^ yr), the CgO and C7O 
abundances are 6.8 x 10^^^ and 1.3 x 10~^° respectively. 
The HCeO and HC7O abundances are 5.8 x 10"^° and 
1.5 X 10~^ respectively. These are comparable to the C 3O 
abundance observed in TMC-1 bv lBrown" eEaDdHH of 
1.4x10"^*^, which suggests that large carbon chain oxides 
may be detectable in TMC-1 and similar carbon-chain- 
rich clouds. The predicted carbon-chain-oxide abun- 
dances depend on the availability of gas-phase 0-atoms 
and are therefore significantly less in the higher-density 
environment of L1544 as a result of depletion of 0-atoms 
onto dust. 
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TABLE 3 

Modeled fractional abundances of 
carbon-chain-oxide neutrals and 

ANIONS 



Model I 



Species 



TMC-1 



L1544 



C30 


2.31 X 10" 


10 


5.77 X IQ- 


ll 


CeO 


6.78 X 10- 


11 


2.35 X IQ- 


12 


C7O 


1.33 X 10" 


10 


1.02 X 10- 


12 


HCeO 


5.81 X 10" 


10 


2.58 X 10- 


12 


HC7O 


1.46 X 10" 


09 


2.16 X 10- 


12 


CeO- 


6.20 X 10" 


12 


1.25 X 10- 


13 


C7O- 


2.08 X 10- 


11 


1.16 X 10- 


13 


HCeO" 


7.38 X 10- 


12 


2.79 X 10- 


14 


HC7O- 


1.45 X 10- 


11 


1.25 X lo- 


14 



Note. — Fractional abundances (with 
respect to n^i ) are from Model III, with 



riHa = 10 cm , 
1, and = lO*' 
L1544. 



lO'' yr for TMC- 
t = 10** yr for 



The most important route to the smaller cumulene ox- 
ide C3O in our model is through the dissociative recombi- 
nation of H3C3O+, which is formed from the reaction of 
C2H;^ with CO. We do not consider the equivalent mech- 
anisms for CeO or C7O synthesis {via C5H3 and CgHj^, 
respectively) , because these ions are much less abundant 
in the model than C2H|^, especially at early times. 

The addition of carbon-chain oxides to our models has 
a relatively minor impact on the abundances of other 
species. They act to boost the abundances of long carbon 
chains and anions by a few percent due to a recycling 
effect, for example: 



CeH- 



followed by 



O 



HCfiO 



(6) 



HCgO + C ^ CeH + CO (7) 

whereas in previous models {e.g. IWalsh et alj [20091 ) . 
0-atom reactions resulted only in the shortening of car- 
bon chains through CO formation: 



CeH- 



O 



C,,H- 



CO 



(8) 



Further dedicated laboratory studies are required in 
order to firmly establish the products of reactions be- 
tween anions and oxygen atoms. In addition, observa- 
tional searches for carbon-chain oxides and their anions 
will help to constrain the relevant branching ratios. The 
neutral (H)C„0 species (for n = 6,7) have known rota- 
tional spectra and possess large dipole moments, making 
them suitable candidates for radio- astronomical searches. 



6. 



TEMPORAL EVOLUTION OF THE 
ANION-TO-NEUTRAL RATIO 



Figure [5] shows the behaviour of the CeH^ anion-to- 
neutral ratios as a function of time at H2 densities of 10*, 
10^ and 10^ cm~^ (red, green and blue curves, respec- 
tively), for each of the three models (i-m)- The anion- 
to-neutral ratio patterns for CeH^ are representative of 
those for the other hydrocarbon anions in the model, 
and show clear variation as a function of cloud age and 
density. The anion-to-neutral ratio curves are distinctly 
different between the three models primarily as a result 




0.001 



Model II 



Model III 



0.0001 




0.001 



1e-05 




1e-05 



Fig. 5. — Solid lines show modelled CgH anion-to-neutral ra- 
tios (left ordinate, linear scale) as a function of time for gas den- 
sities n,H2 = 10* — 10^ cm" ^ in the case of no accretion (Model I; 
top), freeze-out accretion (Model II; middle) and accretion + CR- 
desorption (Model III; bottom). Curves have been truncated when 
neutral CsH fractional abundances fall below 10" Dashed lines 
show the modelled gas-phase atomic oxygen fractional abundances 
(right ordinate, logarithmic scale). 

of their dependence on the gas-phase abundances of elec- 
trons and atomic H, C and O, which themselves vary with 
time and density. A common trend between the three 
models is that in denser media, faster cation recombina- 
tion rates reduce the free electron abundances during the 
early stages of cloud evolution, resulting in lower anion 
production rates. At all modelled densities, the anion- 
to-neutral ratios tend to increase over time as H, C and 
O abundances fall and electron abundances rise. 

For reference, at = 10^ cm~'^ our electron abun- 
dance as a function of time closely matches that shown 
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in Figure 1 of [Roberts et alj (|2003D . In our Model ii 
(Figure m middle panel), the steep rise in the anion-to- 
neutral ratio is mainly a consequence of the rising elec- 
tron abundance that occurs during the advanced stages 
of freeze-out. Dashed lines in Figure [S] show the gas- 
phase 0-atom fractional abundance, which may be used 
as an indicator for the level of depletion. We note that 
our models show qualitative agreement with the trends 
in the CO depletion levels (Sco) observed in TMC-1, 
L1512 and L1544; the CO goes from almost undcpleted 
in TMC-1 to highly-depleted in the denser gas of L1544 
(see Table El). 

Comparing Models i and ii (top and middle pan- 
els of Figure [5]), the effects of gas-grain interactions 
are clear. Gas-phase species (especially lighter, atomic 
species), freeze out onto the dust grains at a faster rate 
in denser media, which has a strong impact on the an- 
ion abundances over time: CrH" react s quickly with 
atomic oxygen (|Eichelberger et aI]|2007D , so the deple- 
tion of O from the gas contributes towards a steady rise 
in [C6H~]/[C6H]. As the depletion of gas-phase species 
proceeds towards total freeze-out - indicated in the fig- 
ures by the rapidly falling 0-atom abundances - electron 
densities rise by a factor of a few. This occurs primarily 
as a result of the falling abundances of HsO"*" and HCO''". 
These species have rapid recombination rates, so their 
loss results in elevated electron densities. Accordingly, 
the [C6H^]/[C6H] ratio rises to ^ 50% towards the end 
of the cloud's lifetime {i.e. towards the time of complete 
freeze-out). 

Such large anion-to-neutral ratios are incompatible 
with the interstellar anion observations made thus far 
(shown in Table [2]). Given the steeply-rising nature 
of the [C6H~]/[C6H] curves in Model ii and the fact 
that the CgH" abundance remains at an observable level 
(> 10"^'^) until relatively late times, a broad range of 
anion-to-neutral ratios from approximately 0-50% is pre- 
dicted. However, all values observed so far have been 
less than 10%. It seems statistically unlikely that obser- 
vations have captured only the youngest clouds on the 
lower parts of these curves. 

Upon inclusion of cosmic-ray-induced desorption of 
atoms in Model ill (Figure [5l bottom panel), the anion- 
to-neutral ratios in clouds with nn^ = 10^ cm~'^ and 
riHs = 10^ cm~'^ obtain values of only a few percent, 
which are consistent with observations. Including des- 
orption in these models keeps the electron density at a 
steady value of ^ lO^^riHa throughout the time of peak 
carbon-chain synthesis until well after the C — > CO con- 
version has occurred. As a result, radiative electron at- 
tachment rates and hence anion abundances remain ap- 
proximately constant (within a factor of a few) for most 
of the cloud's lifetime. At 71H2 = 10"* cm^"^, the results 
match observations at 'early time' (~ 10^ yr), but to- 
wards later times, [C6lI^]/[C6H] continues to rise (due to 
rising electron densities), and levels off at just over 10%. 
This result implies that anion-to-neutral ratios greater 
than presently observed may be detectable in quiescent 
clouds significantly older than TMC-1. The decline of 
carbon-chain abundances at late times will make such 
observations challenging, however. 

7. DISCUSSION 



The CgH^ anion-to-neutral ratio observed in L1527 
(9.3%) is significantly larger than in the other sources 
with similar densities (see Table [2]). Accordingly, Model 
III gives a poorer agreement with L1527 than L1544 and 
L1521F. This discrepancy might be related to the pres- 
ence of the low-mass protostar IRAS 04368-1-2557 at the 
center of L1527, heating from which is not accounted 
for by our model. Fits to the sub-mm SED have shown 
that the dust surrounding IRAS 04368-1-2557 has a tem- 
perature of 22-35 K. The large abundances of carbon- 
chain-bearing molecules (including C„II, C„II~, IIC„N 
and CII3CCII) observed bv'S akai et al.l (|2008 ) have been 
interpreted as arising in the warm protostellar envelope 
as a result of 'warm carbon-chain chemistry'. It has 
been theorized that in this region, methane ice subli- 
mates from the dust grain surfaces and participates in a 
gas-phase chemistry which results in the production of 
la rge abundances of carbo n-chain-bearing species. 

iHarada &: HerbstI (|2008D modelled the gas-phase chem- 
istry of L1527 using a similar approach to our Model i 
(with 71H2 = 10^ cm~^), but with initial abundances ap- 
propriate for a protostellar envelope at 30 K. Volatile 
molecular ices (including CO and CH4) were theorized 
to have been formed on the dust, then released into the 
gas phase during protostellar warm-up . The large initial 
gas-ph ase methane abundance used by IHarada fc Herbs^ 
(|2008D results in greatly-accelerated synthesis of carbon- 
chain-bearing s pecies as a result of ion-m o lecule chem- 
istry (see e.g. IBrown &: CharnlevI 119911 : IHassel et al.l 
|2008[ ). which explains why their modelled pqlyyne and 
cyanopolyyne abundances (jHarada fc Herbs^ 120081 Fig- 
ure 5) are able to match observations after only a few 
thousand years. Our Models II and III, on the other 
hand, give a satisfactory agreement with this source af- 
ter ~ 10^ yr, and show that warm carbon-chain chem- 
istry is not necessarily required to explain the large car- 
bon chain abundances in this source. Instead, the de- 
pletion of oxygen atoms onto dust may be sufficient to 
permit an efficient carbon chain chemistry to develop. 
The presence of abundant hydrocarbon anions is also 
likely to boost the abundances of carbon-chain-bearing 
species (see , Walsh et al..2009.) . However, warm-up chem- 
istry may explain the elevated L1527 anion-to-neutral 
ratio if it can produce greater electron densities or re- 
duced H and 0-atom abundances; to test these ideas 
would require a more detailed chemical model. Indeed, 
Harada & Herbst (200_8) calculated a CgH^ anion-to- 
neutral ratio of about 50% in L1527 (five times greater 
than observed), which is likely to be at least partly a 
result of their low initial H, C and 0-atom abundances, 
which are are several orders of magnitude less than ours. 

Agreement of our model with observations at t ^ 
10* yr is in apparent contradiction with the age of L1527, 
which is likely to be > 10^ yr based on the estimated 
lengt h of the pre-stellar core phase ('^10^ vr: lKirk et al.1 
|2005[ ). However, the gas in the protostellar envelope has 
undergone collapse from lower density to arrive at the 
present value of njja = 10^ cm~'^. Given the density- 
dependence of the timescale for formation of peak carbon 
chain abundances (see Figures [T][3]) , our derived time-of- 
best-agreement cannot be treated as a reliable age indi- 
cator. 

The adopted initial abundances for our models (Table 
[T]) assume the cloud evolves from a neutral state at t = 0. 
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This is contrary to iWalsh et alj (|2009[ ) who assumed all 
the C and S atoms to be initially ionized. As can be 
seen by comparing the results of our riHa = cm~'^ 
model with theirs, the choice of initial ionisation condi- 
tions has a profound impact on the neutral and anionic 
polyyne abundances at very early times {t < 500 yr). 
The large abundance of free electrons in the 'initially- 
ionized' model results in rapid anion production. Conse- 
quently, given the propensity of hydrocarbon anions to 
stimulate synthesis of enhanced polyyne and polyacety- 
lene abundances via anion-neutral reactions, the carbon- 
chain-bcaring species can reach abundances orders of 
magnitude greater than in our 'initially-neutral' mod- 
els. The density-dependence of the chemical timescales 
(see Section [3]) determines whether this 'very early-time' 
peak in carbon chain abundances has any lasting impact 
on the abundances of the species of interest. 

For instance, at tt-Hj = 10^ cm"'^, by the time of best 
fit with observations {t ~ 10^ yr), the abundances in the 
initially-neutral and initially-ionized models have con- 
verged so as to be practically indistinguishable. The 
same can also be said for the models at nn^ = 10^ cm~^. 
However, towards higher densities (n > 10^ cm~'^), con- 
vergence does not occur sufficiently quickly, and the ele- 
vated carbon chain abundances persist through until the 
time of freeze-out. These effects are more pronounced for 
longer carbon chains, and as a result, the CeH abundance 
is about a factor of two greater for = 10® cm"'^ at 
around the time of best fit with observations {t ~ 10^ yr), 
when carbon is initially ionized. Anion-to-neutral ratios 
are also greater by a factor of about 0.5. Although stiU 
in good agreement with observations, it is deemed that 
models with C (and S) initially ionized are less repre- 
sentative of the conditions present in such dense clouds, 
which are generally found in regions of high extinction 
where ionizing UV photons cannot penetrate and the gas 
is therefore neutral. 

Regarding the possible use of anions as chemical 
probes, the reaction rates of CgH" w ith O and H atoms 
were measured in the laboratory by lEichelberger et all 
()200iD as 5.4 x lO"^" and 5.0 x 10-^° cm3s-\ respec- 
tively. Because of this similarity, the dominant CeH" 
destruction mechanisms can be considered to be by re- 
action with H atoms for gas-phase [0]/[H] < 1 and by 
reaction with 0-atoms for [0]/[H] > 1. Although we use 
a slightly larger rate for reactions of CgH" with C-atoms 
(1.0 X 10~^ cm^ s~^), the C-atom abundances remain suf- 
ficiently low in our models that this reaction is never a 
competitive cause of CgH" destruction. The N-atom de- 
struction reactions are also a couple of orders of mag- 
nitude less rapid. At a density n^i^ — 10^ cm"'^, H is 
more abundant than O by a factor of at least 2.5 for all 
times in each of our models. However, at higher densi- 
ties (riHa ^ 10^ cm~'^) - as a consequence of the more 
rapid gas-phase chemistry and H2 formation on dust - 
H becomes up to 2 orders of magnitude less abundant 
than O at early times. In these circumstances, at around 
the times of best agreement of the models with observa- 
tions, the anion abundances are closely dependent on the 
0-atom abundances. Anion abundances are also sensi- 
tive to the number density of free electrons. Therefore, if 
the electron abundance is known, anion-to-neutral ratios 
may be used as a probe of the gas-phase oxygen abun- 



dance in dense interstellar clouds. 

The absolute and relative carbon chain and anion 
abundances in our models are dependent on the gas ac- 
cretion rates, which scale with the total dust grain sur- 
face area and the sticking coefficients, both of which are 
subject to uncertainty. To address this, we have run 
models employing a factor of four range in dust surface 
area (by taking Xq = 2Xc and X'^ = Xg/2), and find 
that although the modelled abundances and times of best 
agreement with observations change slightly (greater dif- 
ferences occurring in the higher-density models), the 
overall results of our study are unaffected. Increased ac- 
cretion rates {X'q = 2Xg) lead to a stronger, earlier 
freeze-out peak in the polyyne and anion abundances, 
combined with a steeper rise in the electron abundances 
and anion-to-neutral ratios at late times. Reduced ac- 
cretion rates {Xq = Xg/2), on the other hand, result in 
significantly greater late-time abundances, allowing com- 
plex molecules to exist in the gas for almost twice as long. 
Examination of the effects of species-dependent differen- 
tial accretion rates, further to what may be inferred from 
our treatment of atomic desorption in Model ill, is be- 
yond the scope of the present study. 

8. SUMMARY 

Similar to prev i ous models for a nions in TMC-1 
()Millar et al.l [20071: IWalsh et all [20091) . for a density of 
riHa = 10'* cm~^ we obtain good agreement between 
modelled and observed anion-to-neutral ratios and ab- 
solute polyyne abundances at times ~ 10^ yr. Simi- 
larly good agreement is found for the quiescent cloud 
Lupus lA. This result is independent of the inclusion of 
gas-grain interactions, which are negligible at this time 
and for this relatively low gas density. Models give a 
poor fit with observations at later times due to their 
falling cyanopolyyne abundances, although the freeze- 
out of atomic oxygen does permit the polyynes to remain 
abundant up to ~ 10® yr. 

The chemical timescale for the formation of abundant 
carbon-chain-bearing species (including the polyynes, 
cyanopolyynes and hydrocarbon anions) is inversely de- 
pendent on the cloud density. Thus, for clouds of mod- 
erate density (riHa — 10^ cm~'^), such as L1512 and 
L1251A, the best fits between modelled and observed 
abundances of carbon-chain-bearing species occur at ear- 
lier times ('-^ 7 X 10'* yr). At this higher density. Models i 
(no accretion) and ill (accretion + CR-desorption) match 
the observations better than Model 11 (accretion with no 
desorption) . 

At densities riHa — 10® cm"'^, the best fit with obser- 
vations of L1544, L1521F and L1527 is for Model ill, 
at t ~ 10^ yr. We find that by accounting for gas- 
grain interactions (in particular, the freeze-out of atomic 
oxygen), the large observed abundances of carbon-chain- 
bearing species in these environments can be reproduced. 
However, a sufficient number of small 0-bearing ions 
are required to remain in the gas-phase to maintain a 
steady electron abundance and therefore keep the radia- 
tive electron attachment rates down, in order to repro- 
duce the relatively low observed anion-to-neutral ratios 
(on the order of a few percent) in LI 544 and LI 52 IF. 
Without atomic desorption, rising electron abundances 
around the time of freeze-out result in excessively large 
anion-to-neutral ratios. 
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In summary, the gas-grain interaction, through its 
influence on the chemistry (in particular, the O and 
H-atom abundances and the electron density), is found 
to have profound effects on the abundances of carbon- 
chain-bearing species and their anions in dense interstel- 
lar clouds with ^ 10^. Due to their destructive reac- 
tions with atomic oxygen, neutral and anionic polyynes 
are predicted to be most abundant in oxygen-depleted 
environments. In addition, 0-atom -|- anion reactions 
are found to have the potential to produce as-yet unob- 
served interstellar species including carbon-chain oxides 
and C2H-. 

The sensitivity of the interstellar CeH" anion-to- 



neutral ratio to the gas-phase O, H and abundances 
could allow anions to be used as a probe of these funda- 
mental chemical species. 
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